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ABSTRACT 

We present a method for determining mean light-weighted ages and abun- 
dances of Fe, Mg, C, N, and Ca, from medium resolution spectroscopy of un- 
resolved stellar populations. The method, pioneered by Schiavon (2007), is im- 
plemented in a publicly available code called EZ_Ages. The method and error 
estimation are described, and the results tested for accuracy and consistency, 
by application to integrated spectra of well-known Galactic globular and open 
clusters. Ages and abundances from integrated light analysis agree with studies 
of resolved stars to within iO.ldex for most clusters, and to within ±0.2 dex for 
nearly all cases. The results are robust to the choice of Lick indices used in the fit- 
ting to within ±0.1 dex, except for a few systematic deviations which are clearly 
categorized. The realism of our error estimates is checked through comparison 
with detailed Monte Ca rlo simulations . Fina lly, we apply EZ_Ages to the sample 



of galaxies presented in iThomas et al.l (120051 ) and compare our derived values of 



age, [Fe/H], and [a/Fe] to their analysis. We find that [ct/Fe] is very consistent 
between the two analyses, that ages are consistent for old {Age > 10 Gyr) popu- 
lations, but show modest systematic differences at younger ages, and that [Fe/H] 
is fairly consistent, with small systematic differences related to the age system- 
atics. Overall, EZ_Ages provides accurate estimates of fundamental parameters 
from medium resolution spectra of unresolved stellar populations in the old and 
intermediate-age regime, for the first time allowing quantitative estimates of the 
abundances of C, N, and Ca in these unresolved systems. 

Subject headings: methods: data analysis, galaxies: abundances, galaxies: star 
clusters 
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Introduction 



The study of stellar populations in galaxies is going through very exciting transforma- 
tions. The advent of extensive and high-quality data sets for galaxies in the local and distant 
universe poses a demand for ever improving stellar population synthesis models. In this pa- 
per, we present a new tool, named EZ_Ages, that simultaneously address two on-going needs 
in the field of stellar population synthesis: 1) it extracts accurate ages and abundances of 
several elements from integrated spectra of stellar populations, and 2) it performs this task 
in an automatic fashion, which is well-suited for applications to large data sets. 

Early stellar population synthesis models for single-burst populations — so called "single 
stellar populat i ons" or SS Ps — were based on lo w resolution stellar spectra. The models of 
Worthey et al.l (119941 ) and IVazdekis et al.l (Il996l ) used optical absorption- l ine in dices with 8- 
10 A resolution and limited spectral coverage, while ISruzual Sz CharlotI ( 1993 ) synthesized 
model spectra with broad wavelength coverage (ultraviolet through infrared) and resolution 
of 10-20 A in the optical. All models were for solar- abundance ratios onl y. Since then, a new 
gener atio n of SSP models have pushed to higher resolution (2-3 A for 



Weiss et al. 



Maraston 



1995 



2OO5I : 



Tantalo et a 



2OO3I and IVakdekis et al.l 12003) and to variable abundance patterns (iBorges et al 



Schiavon 2007 and Coelho et al. 



19981: iTrager et a. 



2003). 



Bruzual fc Charlo 



2000l : lThomas. Maraston. fc Bendei]l2003 : 



1995 



It is well known that many stellar populations, including the Galactic halo and bulge. 
Galactic globular clusters, massive early-type galaxies, and the bulges of spiral galaxies, have 
super-solar [a/Fe]. In resolved systems, where abundance analysis of individual stars is possi- 
ble, it has been known for a while that not all a elements are equally enhanced. On the other 
hand, in unresolved systems it has been only recently that evidence started accumulating 
to the fact that observed stellar populations have abundance patterns which cannot be de- 
scribed by only two parameters (ie. [Fe/H] and [a/Fe]) but rather require more complicated 
abundance patterns. In particular , it seems that not all the a elements are equally en- 



hanced in early-type ga l axies (e.g.. IVazdekis et al. 



1997 



Henry fc Wortheyl 1 19991 : iThomas. Maraston. fc Bender 



Worthev 



2003 



1998 



Trager et al.lll998 



SchiavonI 120071 ). In the best 



studied, closest, resolved spheroidal system, the Gal actic bulge, there is also recent evi- 
dence that not all g -elements are equally enhanced (IFulbright. Rich. &: McWillianu 12005 



Fulbright et al.ll2007l ). In particular, IFulbright et al.l (|2007l ) found that, while Mg is strongly 



enhanced relative to Fe ([Mg/Fe] ~ +0.3 across a range of [Fe/H]), other a elements such 
as Si, Ca, Ti, and O tend to follow a disk-like trend, with [X/Fe] decreasing with increas- 
ing [Fe/H]. Interestingly, in the case of oxygen — ar guably the most important among the 
a elements — they found [0/Fe] ~ at [Fe/H] = 0. iBensby. Feltzing. fc LundstromI (l2004r) 



show a similar result for Milky Way disk stars in both the thin and thick disk. lAthey 
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( 120031 ) demonstrates that warm ionized gas in a sample of 7 local early-type galaxies has 
mean [0/H] ^ —0.1, while the stellar populations of the galaxies have super-solar [Mg/H]. 
Clues to variations of the abundances of other elements, such as C and N, particularly 
as a function of environment, have been found in a number of studies. The Lick indices 
C24668 and CN2, which are sensitive to C and O (and N, in the case of CN2), have been 
shown to be weaker in early-type galaxies in dense environments as compared to galax- 
ies in t he field, while the strength o f the Mg b index remains unchanged across environ- 



ments ( ISanchez-Blazquez et al.l 120031 ) . Finally, a number of studies have found evidence 



that Ca seems to behave as an Fe-£ 



ing the trend of 


;her a elements ( 


Vazdekis et al. 


199' 


'': Worthev 


1998: 


Traeer et al. 


1998: 


Henrv & Wort 


ley 


1999: 


Thomas. Maraston. & Bender 


2003). Thi 


3 result has recently been 
>how that the Ca4227 in- 


questioned by 


Prochaska et al. 


( 


2005) and 


Schiavo] 


Q (2 


DOtI), who £ 



dex is strongly affected by CN lines, which may be partially responsible for the low [Ca/Fe] 
suggested by other authors. Although a-enhanced models are an important improvement 
over solar-scaled-only, it is fundamental that models with more degrees of freedom be de- 
veloped, which are able to explore the complexity suggested by the data. By estabhshing 
a more detailed snapshot of the abundance patterns of stars in galaxies, these models can 
pose fundamental constraints on models of galaxy formation. 



SchiavonI (120071 ) presents a new set of SSP models which do exactly that. These models 
explicitly include adjustable abundance patterns for multiple elements, allowing the user to 
separately vary the abundances of C, N, O, Mg, Ca, Na, Si, Cr, and Ti. This paper describes 
an algorithm for searching through the space of variable abundance models to find the 
abundance mix which best fits the Lick indices measured in a given input stellar population 
spectrum. This algorithm has been implemented in an IDL code package, called "EZ_Ages" 
which is publicly available for use. Section [2] describes the sequential grid inversion algorithm. 
In ^ we present ages and abundances estimated by EZ_Ages for the Galactic globular 
clusters NGC 6121, 47 Tuc, NGC 6441, and NGC 6528, as well as the open cluster M67. 
These values are compared to results in the literature determined from isochrone fitting of 
cluster photometry and from high-resolution spectral abundance analysis for individual stars 
in the clusters. In §1] we apply EZ_Ages to galaxy data from the literature, contrasting our 
results with those by other authors. Section [5] summarizes our conclusions. 

The beta-version of EZ_Ages is available for download, along with instructions on in- 
stallation and use 111 



1 www . ucolick. org / '^graves /EZ_Ages . html 
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2. Sequential Grid Inversion Algorithm 
2.1. Schiavon Simple Stellar Population Models 



m 



The code EZ_ Ages is designed for use with the single stellar population models described 



SchiavonI (120071 . hereafter S07). The details of the model are described in SOT; here we 



briefly discuss the aspects of the model relevant to its practical use. 

These models include a ch oice of two se t s of i sochrones by the Padova group: one 
is the solar-scale d isochrones of iGirardi et al.l (120001 ) and the other is the a-enhanced set 
of isochrones by ISalasnich et al.l (120001 ) with average [a/Fe] = +0.42 and with individual 
abundance ratios chosen to match metal-poor field stars^l However, th e mod els based on 
the a-enhanced isochrone must be used with caution because Iweiss et al. ( 2006 ) have shown 
that, due to a problem with the opacity tables adopted, those isochrones predict temperatures 
that are slightly too high at both the main-sequence turnoff and the giant branch. The effect 
is strongest at solar metallicity, where it results in ages that are slightly too old. For details, 
see §4.3 of S07. Developing isochr ones with variable abundance ratios is an active area of 
current research (see recent work by lDotter et al.ll2007l ) and we hope to incorporate a broader 
range of isochrones into EZ_Ages in the near future. 

In addition to a choice of isochrones, the S07 models allow the user to select an input 
abundance pattern. The abundance ratios that can be varied are: [Mg/Fe], [C/Fe], [N/Fe], 
[0/Fe], [Ca/Fe], [Na/Fe], [Si/Fe], [Cr/Fe], and [Ti/Fe]. The S07 models output predicted 
Lick index measurements for a grid of age and [Fe/H] values, based upon the input abundance 
ratios. The models span the age range 0.1 < t < 15.8 Gyr, with —1.3 < [Fe/H] < +0.2 for 
the solar-scaled isochrone and —0.8 < [Fe/H] < +0.3 for the a-enhanced isochrone, though 
only models older than 0.8 Gyr are computed with variable abundance ratios. 

The Lick indices included in the S07 model are B.6a, HSp, H'Ja, Hif, H/?, CNi, CN2, 
Ca4 227, G4300, Fe4383, C .^4668, Fe5015, Mg?, Mg 6, Fe5270, and Fe533 5, as defined in Table 
1 of lWorthev et all Jl994h and Table 1 of IWorthev fc Ottavianil Jl997h . The output model 
indices can then be compared with Lick indices measured in the integrated spectrum of a 
real stellar population (i.e., a stellar cluster or a galaxy) to determine the best-fitting mean 
light-weighted age and [Fe/H] for the stellar population, for a given set of model abundance 
ratios. 



^The following a elements are enhanced as shown: [0/Fe] = +0.50, [Ne/Fe] = +0.29, [Mg/Fe] +0.40, 
[Si/Fc] = +0.30, [S/Fe] = +0.33, [Ca/Fe] = +0.50, [Ti/Fe] = +0.63, [Ni/Fe] = +0.02. All other abundance- 
ratios are solar. 
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The Lick indices included in the S07 model are summarized in Table [H along with their 



main abun dance sensitivities. These are taken from the works of IServen. Worthev. fc Brilev 



( 2005 ) and lKorn. Maraston. &: Thomas! (120051 ) . The sensitivities from lKorn. Maraston. fc Thomas 



( 120051 ) are taken from their Tables 4 and 5, which give abundance sensitivities for turnoff stars 
and giant branch stars respectively. We have chosen to include only those two evolutionary 
phases in the table, as they dominate the stellar population model spectrum (accounting 
for ~ 90% of the integrated light), even though all three evolutionary phases are taken into 
account in the model computations. In addition to the sensitivities to individual element 
abundances shown in the table, all indices are sensitive to changes in the total metallicity. 
We interpret the effect of [Z/H] on Fe5270 and Fe5335 as being due primarily to changes in 
Fe because none of the other elements investigated by KMT05 appear to affect these indices. 
Table [T] shows reasonable agreement between the two different sets of sensitivity tables. 

Abundance variations enter the model in two ways, firstly through the choice of isochrone 
(solar-scaled or a-enhanced). At this stage, the model takes into account the effect of the 
abundance mix on stellar evolution — specifically, how the chemical mix affects temperatures 
and luminosities of stars of different masses and at different evolutionary stages. Admittedly, 
this is done in an approximate fashion, given that unfortunately we only have isochrones com- 
puted for two abundance patterns. Secondly, the effect of individual elemental abundances 
on line and continuum opacities in stellar atmospheres i s incorporated into model Lick inde x 



predictions, using the response functions computed by iKorn. Maraston. fc Thomas! (120051 ). 
At this second stage, it is possible to vary elemental abundances independently; at the first 
stage, the only choices are [a/Fe] = and +0.42, with total metallicity varying. 

It should be stated clearly that these models are cast in terms of [Fe/H], and not 
total metallicity, [Z/H]. This reflects our choice to deal explicitly with quantities that can be 
inferred from measurements taken in the integrated spectra of galaxies — total metallicity not 
being one of them, given our current inability to use integrated spectra of stellar populations 
to constrain the most abundant of all metals, oxygen (see discussion in §4.4.1 of SOT). 
Another advantage with casting models in terms of [Fe/H] is that each elemental abundance 
can be treated separately, so that the effect of its variation can be studied in isolation from 
every other elemental abundance (at the cost of varying the total metallicity). In the case 
of models cast in terms of [Z/H], it is impossible to vary the abundance of a single element, 
because enhancing one element means decreasing the abundances of all other elements to keep 
total metallicity constant. As an example, in models cast in terms of [Z/H], a solar-metallicity 
([Z/H] = 0) a-enhanced population has only a slightly higher abundance of a elements than 
the Sun, while having a much lower iron-abundance. It is the large iron abundance difference 
that is responsible for the bulk of the difference between those two hypothetical models, not 
the small difference in a-element abundances. For this reason, claims that higher order 
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Balmer lines are strongly affected by a-element abundances ( iThomas et al.l 120041 ) should be 



taken with caution (see §4.3 SOT for a thorough discussion of this point). 



2.2. Fiducial Age and [Fe/H] 

In practice, one would like to use absorption line measurements in the integrated spec- 
trum of a given stellar population to determine not only the age and [Fe/H] for a fixed set 
of abundance ratios, but to find the best set of input abundances as well. The brute-force 
method of finding the best-fitting abundances would involve generating a set of models that 
span the available space of age, [Fe/H], and abundance patterns, and then identifying the 
model whose predicted Lick indices best match those measured in the spectrum. However, 
with nine variable abundance ratios, the search space quickly becomes very large. If each of 
those abundance ratios is represented by only four possible values, with 23 model ages and 
four possible values of [Fe/H], that would mean creating more than 24 million models! 

Instead, we choose to perform a directed search for the best model, taking advantage 
of the fact that various Lick indices are sensitive to only a few, often different, elemental 
abundances, as indicated in Table [TJ These sensitivity variations have been modeled by 
several gr oups using spectrum synthesis, ba sed on model stellar atmospheres ( see recent 



papers by iKorn. Maraston. fc Thomas! |2005| and IServen. Worthey. &: Brileyl |2005| ) . A more 



detailed discussion of the motivation behind this method can be found in SOT. 

EZ_Ages begins by computing a set of models with the chosen isochrone (either solar- 
scaled or a-enhanced), using solar abundance ratios for the stellar atmospheres. It then uses 
a pair of lines that are sensitive to age and [Fe/H], but relatively insensitive to other elemental 
abundances — a Balmer line and an Fe line — to determine a fiducial age and [Fe/H] for the 
stellar population in question. From Table 1 we can see that H/5 and the iron lines Fe52T0 
and Fe5335 are good choices for determining the fiducial, as they are mostly insensitive to 
other element abundances. In EZ_Ages, the default choice of indices for calculating fiducial 
age and [Fe/H] are B.j3 and (Fe), an average of Fe52T0 and Fe5335. 

The top left panel of Figure [1] shows a plot of model grids for H/3 and (Fe). Dotted lines 
connect constant-age models (from top to bottom: 1.2, 2.2, 3.5, T.O, and 14.1 Gyr) while solid 
lines connect models with the same [Fe/H] (from left to right: -1.3, -O.T, -0.4, 0.0, and +0.2). 
The square shows a sample data point for a galaxy with age ~ T Gyrs and [Fe/H] ^ —0.2. 
In this figure, it is easy to see which box of the grid encloses the data point in question. 
This gives a bounded range in age and [Fe/H] for the data point. A two-dimensional linear 
interpolation within the gridbox is then used to convert the data point in index-index space 
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into fiducial values for the age and [Fe/H]. 

2.3. Determining Non-Solar Abundance Ratios 

Once a fiducial age and [Fe/H] have been determined from the H/5-(Fe) model grid, one 
can then make a similar grid plot with an index which is sensitive to non-solar abundance 
ratios, for instance a Balmer line plotted against Mg b. Mg b is strongly affected by the 
[Mg/Fe] abundance ratio, at fixed [Fe/H]. If the chosen S07 model is a good match to the 
data, then the age and [Fe/H] estimated by grid inversion from the H/5-Mg b plot should 
match the fiducial age and [Fe/H] from the H/3-(Fe) plot. In our example in Figured], this is 
clearly not the case; the H/3-Mg b plot in the top right panel of Figure [1] gives a larger value 
of [Fe/H] than the H/3-(Fe) plot in the top left panel. This is an indication that the example 
galaxy has a larger [Mg/Fe] ratio than is predicted in a solar-scaled model. The sequential 
grid inversion algorithm then increases the input value of [Mg/Fe] and recomputes the S07 
model. Increasing [Mg/Fe] will increase the strength of the Mg b index at fixed [Fe/H], which 
will have the effect of sliding the entire set of model grids to the right in the H/5-Mg b plot, 
while keeping them essentially unchanged in the H/3-(Fe) plane. This lowers the value of 
[Fe/H] that is estimated by grid inversion in the H/5-Mg b plot, bringing it into agreement 
with that obtained in the H/?-(Fe) diagram. This can been seen in the lower panels of Figure 
m which show H/3-(Fe) and H/9-Mg b plots for an S07 model computed with [Mg/Fe]= +0.3. 
For the data shown in this plot, the age and [Fe/H] values indicated by the model grids are 
consistent for both metal lines. In the sequential grid- inversion algorithm implemented in 
EZ_Ages, [Mg/Fe] is increased incrementally until the [Fe/H] estimated from the H/5-Mg b 
plot matches the fiducial [Fe/H] from the H/5-(Fe) plot. This process is then repeated for 
other Lick indices and other abundance ratios. 

The key to the sequential grid inversion algorithm is to proceed with the abundance 
fitting in such a way as to only adjust one abundance at a time. Once fiducial values for the 
age and [Fe/H] of the system have been determined, the goal is to next fit a Lick index that 
introduces only one additional abundance. Mg 6 is a good choice because it is dominated by 
Mg and Fe (see Table 1), with the latter being determined from the H/3-(Fe) grid inversion. 
By adjusting [Mg/Fe] until the H/3-Mg b grid inversion matches the fiducial values, EZ_Ages 
can find the best-fitting Mg abundance. 

Another index which EZ_Ages can use to fit a single element is C24668, which is domi- 
nated by C and O. Oxygen is notoriously difficult to estimate in old (ie. non-star forming) 
unresolved stellar populations, therefore EZ_Ages does not try to model [0/Fe] (but see 
discussion in §3.21) . Instead, it allows the user to choose an input oxygen abundance. For 
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consistency, this should be chosen to match the input isochrone: [0/Fe] = 0.0 for the solar- 
scaled isochrone, [0/Fe] = +0.5 for the a-enhanced isochrone. With [0/Fe] set by the user, 
this leaves only C as a variable abundance which contibutes strongly to C24668. EZ_Ages 
adjusts [C/Fe] so that the [Fe/H] estimated in the H/9-C24668 plot matches the fiducial. 

The only indices which are strongly affected by nitrogen and which might therefore 
allow a fit of [N/Fe] are the two CN lines, CNi and CN2, both of which respond strongly to 
changes in C, N, and O. As discussed above, [0/Fe] is set by the user. This leaves C and 
N as variable abundances contributing to the CN indices. However, if [C/Fe] has already 
been determined by fitting C24668, then the CN indices can be used to fit for [N/Fe] at the 
calculated value of [C/Fe]. The computed value of [N/Fe] is therefore dependent upon the 
calculated value of [C/Fe] and a robust error analysis should propagate errors from [C/Fe] 
into the the errors in [N/Fe]. 

The S07 models include one Lick index dominated by Ca and can therefore be used to 
fit [Ca/Fe]: Ca4227. Because the blue continuum region used in measuring Ca4227 overlaps 
with a CN absorption band, this index is also very sensitive to [C/Fe] and [N/Fe], so that 
these abundance ratios must be fit properly before using Ca4227 to determine [Ca/Fe]. As 
with [N/Fe], errors in the calculated values of [C/Fe] and [N/Fe] contribute to the error in 
[Ca/Fe]. Again, this should be taken into consideration in a complete error analysis (see 

ED- 

The set of available Lick indices and their element abundance sensitivities thus prescribe 
an order for abundance fitting: [C/Fe], then [N/Fe], then [Ca/Fe]. Note that Mg b, the 
preferred index used to fit [Mg/Fe], does not depend strongly on these abundances, nor do 
the other abundances depend on it, therefore it can be fit in any order. 

Based on Table [H it is clear that some indices are preferable for fitting: those that 
are dominated by the elemental abundances that we are using them to estimate and no 
others. Because of this, EZ_Ages has a built-in hierarchy which it uses to choose the indices 
used in abundance fitting. Among the Balmer lines, the preferred index is H/3 because it 
is mildly sensitive to total metallicity, but relatively insensitive to individual abundances. 
It is followed by HSp and B.jp, which are somewhat sensitive to individual elements, and 
lastly H^^ and Hja, which are broader indices (designed to measure Balmer strengths in 
populations with a substantial A star component) and therefore include a larger number of 
other absorption features in the index and pseudo-continuum bandpasses. 

For iron indices, the order of preference is (Fe) (an average of Fe5270 and Fe5335), 
Fe5270, Fe5335, Fe4383, and Fe5015. The indices Fe5270 and Fe5335 are dominated by Fe 
and total metallicity effects, making them relatively clean indicators of [Fe/H]. Of the other 
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two lines, both of which are influenced by element abundance ratios, we prefer Fe4383 over 
Fe5015 for two reasons: Fe5015 falls on a CCD defect in the globular cluster data used in 
^to test the output of EZ_Ages and thus can only be tested for M 67, where it appears to 
perform similarly well to Fe4383, and the Fe5015 bandpass overlaps the [Olll]A5007 emission 
line from ionized gas sometimes present in galaxies, thus limiting its usefulness in stellar 
population modeling of galaxies. 

For fitting [Mg/Fe], the preferred index is Mg b followed by Mg2, which is more sensitive 
to C than is Mg b. For [C/Fe], it is C24668 followed by G4300, for reasons discussed in detail 
in §3.21 For [N/Fe], it is CNi followed by CN2, although this is fairly arbitrary as there is no 
strong reason to prefer one CN index over the other. For [Ca/Fe], the only available index 
is Ca4227. 

It is possible to manually override this order of preference if a given index is unavailable, 
or to investigate the effects of using other lines in the fit. Note that EZ_Ages will only fit an 
abundance ratio if at least one of the indices needed to estimate it is available. Also, abun- 
dances which depend on other abundance ratios will not be fit if the prerequisite abundance 
is not available (e.g., EZ_Ages will not use a CN index to fit [N/Fe] unless [C/Fe] has been 
determined from C24668 or G4300). 

Thus far, we have focused only on the dominant abundances affecting each index, how- 
ever many indices have a weak dependence on other abundances. Particularly if the most 
desirable Lick indices are not available due to limited wavelength coverage, bad skyline sub- 
traction, or nebular emission within the galaxy and other lines have to be used (such as B.6f 
in place of H/? or Fe4383 in place of Fe5270 or Fe5335), the computed fiducial values for age 
and [Fe/H] or other abundance values may change once the abundance fitting is complete, 
because some of the newly determined abundances may affect the indices used to determine 
those fiducial values. To ensure that the final abundance ratio is consistent for all of the 
Lick indices used in the fit, EZ_Ages iterates the fitting process up to four times. It uses the 
best-fitting abundance pattern from the previous iteration to compute new fiducial values 
for age and [Fe/H], then repeats the fitting process until further iterations do not improve 
the fit. 

The S07 models account for variations in Lick index strength due to abundance varia- 
tions in the elements Fe, Mg, C, N, Ca, O, Na, Si, Cr, and Ti. However, only the first five of 
these elements can be fit using the algorithm described above. The elements O, Na, Si, Cr, 
and Ti do not dominate any of the Lick indices modeled by S07 and therefore cannot be fit in 
this process. The EZ_Ages code allows the user to supply a value for any of these unmodeled 
elements to use as input to the S07 models. In the absence of user-supplied values, EZ_Ages 
has a set of default abundances for these elements. By default, it sets [0/Fe] to match 
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the chosen isochrone ([0/Fe] = 0.0 for the solar-scaled isochrones, [0/Fe] = +0.5 for the 
a-enhanced isochrone). The default for the a-elements Na, Si, and Ti is that they are set to 
follow Mg, so that [Na/Fe] = [Si/Fe] = [Ti/Fe] = [Mg/Fe]. The iron-peak element Cr is set to 
follow Fe, so that [Cr/Fe] = 0.0. All of the modeling reported in this analysis was done using 
the default settings for the unmodeled elements and with the solar-scaled isochrone. To see 
the effect of using the a-enhanced is ochrone or super-sol ar [0/Fe] on abundance fitting, see 



If errors in the Lick index data are provided as input, EZ_Ages uses them to estimate 
the uncertainties in the ages, [Fe/H], and abundance ratios calculated by the sequential grid 
inversion algorithm. Uncertainties are assumed to be dominated by measurement errors in 
the line strengths. Systematic uncertainties in the models are ignored, except that the usual 
caveat applies: comparisons of relative ages and abundances between multiple objects are 
more reliable than absolute age or abundance estimates. 

Uncertainties in age and [Fe/H] due to line strength measurement errors are computed by 
displacing the measured data point by the input index error, then redoing the grid inversion 
on the displaced data point to determine the range in age or [Fe/H] constrained by Icr error in 
the measured Balmer or Fe line strength. Because the grid spacing is non-linear (particularly 
in age), errors in the positive and negative direction are computed separately. It has been 
noted by many previous authors that the errors in SSP age and [Fe/H] as inferred from index- 
index plots are not independent due to the fact that the model grids are not orthogonal. As 
can be seen from an inspection of Figure [H age and [Fe/H] errors are correlated in the sense 
that lower inferred ages will lead to higher inferred [Fe/H], and vice versa. EZ_Ages does not 
explicitly model this aspect of the inferred errors in age and [Fe/H] values, but users should 
keep this effect in mind when interpreting results from SSP index-index fitting. 

Any error in the choice of fiducial age and [Fe/H] will result in errors in the abundance 
ratios estimated by fitting to incorrect values of the age or [Fe/H]. Thus uncertainties in 
the fiducial age and [Fe/H] must be propagated through to determine their effect on the 
abundance ratio estimates. 

For example, EZ_Ages assumes that the uncertainty in [Mg/Fe] is dominated by two 
sources: the measurement error of the Mg-sensitive Lick index used in the abundance ra- 
tio determination (e.g., Mg b) and the uncertainty in the fiducial [Fe/H] which EZ_Ages 
attempts to match in the H/3-Mg b plot. In the following discussion, we will use AMg b 




2.4. 



Error Calculations 
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and A[Fe/H] to represent the measurement error in Mg b and the uncertainty in the fiducial 
[Fe/H], respectively. The uncertainty in [Mg/Fe] due to Mg b measurement error and fiducial 
uncertainty are denoted in turn by A[Mg/Fe]Mgb and A [Mg/Fe] gd- 

To determine A[Mg/Fe]Mgb, EZ_Ages recomputes the [Mg/Fe] abundance for Mg &+AMg b 
and Mg &— AMg b, using the larger of these two uncertainties as the estimated A[Mg/Fe]Mgb- 

To determine A[Mg/Fe]fid, EZ_Ages recomputes [Mg/Fe] using the measured vahic of Mg b 
to fit [Fe/H] + A [Fe/H] and [Fe/H]— A [Fe/H], again choosing the larger uncertainty for es- 
timating A[Mg/Fe]fid. These two sources of error are treated as independent and the total 
A[Mg/Fe] is determined by summing in quadrature A[Mg/Fe]Mgb and A[Mg/Fe]fid. 

The process for computing A[C/Fe] is identical to that for [Mg/Fe], substituting a C- 
sensitive Lick index for Mg b. Because [N/Fe] is determined using a CN-sensitive index 
such as CNi, A[N/Fe] depends not only on the CNi line strength measurement error and 
A[Fe/H], but on A[C/Fe] as weU. The contributions to A[N/Fe] due to ACNi and A[Fe/H] 
(i.e., A[N/Fe]cNi and A[N/Fe]fid) are computed as above. An additional contribution due to 
A[C/Fc] is determined by matching the measured CNi line strength to the fiducial [Fe/H], 
using models with carbon abundances given by [C/Fe]-|-A[C/Fe] and [C/Fe]— A[C/Fe]. These 
three sources of error are summed in quadrature to obtain A [N/Fe]. A similar process is 
repeated for [Ca/Fc], with the sole difference that errors in [N/Fe] also contribute to the 
[Ca/Fe] error budget. 

To check that this analysis provides a reasonable estimate of the actual statistical errors, 
we have used a Monte Carlo method to generate a large set of artificial "repeat measure- 
ments" of the Lick indices for a test stellar population and compared the resulting distri- 
bution of ages, abundances, and abundance ratios to the errors calculated by EZ_Agcs. For 
our test data, we used Lick index values measured in a stack of ~1 0,000 red sequence galaxy 
spectra taken from the Sloan Digital Sky Survey. As modeled by EZ_Ages, these correspond 
to a stellar population with a mean light-weight ed age of 9.1 Gyr, with [Fe/H]= —0.15, 
[Mg/Fe]= +0.17, [C/Fe]= +0.17, [N/Fe]= +0.07, and [Ca/Fe]= 0.00. 

We ran four Monte Carlo simulations for different error ranges (effectively different S/N 
spectra), in which we assumed the errors were 2%, 5%, 10%, and 20% of the measured index 
values. For each run, we simulated 100 repeat measurements of the data by choosing values 
for each Lick index drawn from Gaussian distributions centered on the "true" data values 
with widths characterised by the assumed 2%, 5%, 10%, and 20% errors. This produced 
a simulation of 100 reahzations for each of the four assumed error ranges. EZ_Ages was 
used to determine the age, [Fe/H], and abundance ratios for each simulated measurement. 
The distribution of these values for the 100 simulated repeat measurements can be compared 
with the error estimate made by EZ_Ages for the original data values and index measurement 
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errors. 

This comparison is shown in Figure [2J The errors reported by EZ_Ages are shown on 
the X-axis, while the y-axis shows the standard deviation a of the distributions of ages and 
abundances produced in each Monte Carlo simulation. As discussed earlier in this section, 
EZ_Ages calculates separately errors in the positive and negative direction. There are thus 
two error values produced by EZ_Ages plotted against each Monte Carlo a. The solid line 
shows a one-to-one relation. The dashed lines outline error estimates that are within ±30% 
of the Icr Monte Carlo spread. 

For the simulations with 2%, 5%, and 10% index measurement errors, the age and 
abundance ratio error estimates from EZ_Ages are all within 30% of the results of the Monte 
Carlo simulation with the exception of the error in the age, which is overestimated by 
about 0.3 Gyr in the 2% error case. Although this is a substantial fraction of the simulated 
H7i? age error, it is well within the uncertainties in the modeling proc ess, which are o f order 



1 Gyr in the ~ 9 Gyr age range corresponding to the test data point (jSchiavonll2007l ). In all 
other cases, the EZ_Ages error estimate is an excellent match to the la errors produced by 
the full Monte Carlo realization. For the simulations with 20% index measurement errors, 
EZ_Ages tends to slightly overestimate the errors as compared to the Monte Carlo simulation, 
but in all cases is still within 35% of the simulated error. 

The good match shown in this test between the error estimates generated by EZ_Ages 
and those produced by a Monte Carlo simulation of the errors suggests that EZ_Ages does in 
fact do a reasonable job of estimating the errors in the age, [Fe/H], and abundance ratio values 
calculated. This includes the propagation of errors between various stages of the modeling 
process. As the data become noisier, EZ_Ages has a tendency to slightly overestimate the 
actual uncertainty in the age, [Fe/H], and abundance due to measurement errors. Figure [2] 
also demonstrates that in order to estimate ages, [Fe/H], and abundance ratios within ±0.1 
dex for a typical early type galaxy. Lick indices must be me asured with ±5% accuracy, which 



requires spectra with S/N ~ 100 A'^ fICardiel et al.lllQQSh . 



2.5. Correlation of Errors in Derived Parameters 



The Monte Carlo error simulation also allows us to investigate the effect of correlated 
errors in the EZ_Ages analysis. Because lines of constant age and [Fe/H] are not orthogonal 
in H/3-(Fe) space, errors in either index measurement will result in anti-c orrelated errors in 
the inferred age and [Fe/H], as discussed in detail in iTrager et al.l (I2OOOI ). For example, an 
error in H/3 that causes the determined SSP age to be too old will cause the measured [Fe/H] 
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value to be too low, as can be seen from the grids in Figure [H Likewise, the sequential grid 
inversion process for fitting abundance ratios depends on the fiducial age and particularly 
on the fiducial [Fe/H] determination. This will tend to produce an anti-correlation between 
[Fe/H] and the measured abundance ratios. 

As shown in Figure [21 the simplified error propagation implemented in EZ_Ages does 
an excellent job of tracking the total error in each stellar population parameter throughout 
the fitting process. However, it does not produce a full error ellipse illustrating correlated 
errors. 

Figure [3] uses the results of the Monte Carlo error simulation to illustrate the correlated 
errors that result from sequential abundance fitting. The 100 simulated data realizations for 
each of the assumed 2%, 5%, 10%, and 20% index measurement errors (yellow, red, blue, 
and green points, respectively) effectively trace out the error ellipses in stellar population 
parameter space. Figure [3^ shows the well-known error anti-correlation between the age and 
[Fe/H] measurements, as expected from the non-orthogonality of the age and [Fe/H] model 
grids in H/3-(Fe) space. 

Panels b-e show the abundance ratio determinations for [Mg/Fe], [C/Fe], [N/Fe], and 
[Ca/Fe] against [Fe/H]. The error ellipses for these data are expected to show an anti- 
correlation of errors because an over-estimate of [Fe/H] will tend to result in an under- 
estimate of the other abundances. For example, if an error in the (Fe) index causes the 
fiducial [Fe/H] determination to be too high, the measured abundance ratios for other ele- 
ments will tend to be too low (at a higher [Fe/H], a lower value of [Mg/Fe] is needed to match 
the same Mg b index). This anti-correlation is observed between [Fe/H] and the abundances 
[Mg/Fe], [C/Fe], and [Ca/Fe]. Interestingly, there appears to be no correlation in the errors 
in [Fe/H] and [N/Fe]. Because [N/Fe] is fit using a CN index, it depends strongly on the 
value of [C/Fe] previously determined by fitting either C24668 or G4300 and thus we expect 
a anti-correlation of errors between [C/Fe] and [N/Fe]. This is indeed the illustrated 
in Figure [3]F, and may explain the lack of correlation in the errors in [N/Fe] and [Fe/H]. As 
in the case of [N/Fe], the fitted values of [Ca/Fe] depend on other abundance fits and thus 
may produce correlated errors. These are illustrated in Figure [3^,h. 

The effect of correlated errors can produce spurious trends in a dataset, thus the corre- 
lations illustrated here should be kept in mind when interpreting abundance patterns from 
EZ_Ages or a similar abundance fitting process. Of course, the best guard against these 
effects is to perform abundance fitting only on spectra with S/N > 100 A^^, corresponding 
to 5% index errors (red points in Figure [3]) or better. 
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3. Testing the Models on Cluster Spectra 

The validation of stellar population synthesis models depends crucially on their abil- 
ity to match the integrated properties of simple systems, such as stellar clusters, which 
can be reasonably approximated by single stellar populations. In fact, this simple tru- 
ism disguises the considerable complexity of the task, as attested by the vast literature 
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; and references therein). S07 



presented a lengthy discussion of the comparison of his models with integrated data for 
well-known Galactic open and globular clusters. Models were computed for a range of ages 
and metallicities, adopting the known abundance patterns of M 5, 47 Tuc, NGC 6528, and 
M 67. These were then compared with Lick index measurements taken in the (very high 
S/N) cluster spectra. As a result, it was found that predictions for Balmer and metal lines 
agreed with the data to within 1-a index measurement errors. Also very importantly, the 
best-matching ages differed by no more than ~1 Gyr from the value obtained from analysis 
of cluster color- magnitude diagrams (CMDs). The only exception was M 5, for which the 
age found was too low, due to the presence of blue horizontal branch (HB) stars, which are 
not accounted for by the models (see §3.21 for a further discussion of this point). 

In this section we ask the inverse question: how well do the models fare when used 
to convert a given set of index measurements into ages and elemental abundances? Again, 
we rely on data for clusters whose abundances and ages are well constrained from analysis 
of individual stellar spectra and deep cluster CMDs. Two important model checks are 
performed. We start by submitting the models to the most challenging of the two tests 
( §3.ip . which consists of checking whether the results obtained with EZ_Ages match those 
from detailed resolved analysis of high-resolution spectra and deep photometry of stellar 
members. While our goal is to match ages and the abundances of Fe, Mg, C, N, and Ca 
to within 0.1 dex our requirement is that differences be smaller than 0.2 dex. Moreover, 
this set of specifications should be met in a quick and non-interactive fashion, suitable for 
applications to large data sets through use of EZ_Ages. As a second check ( §3.2p . we verify 
the internal consistency of EZ_Ages, by comparing the results obtained for a given set of 
observations, when using different combinations of Lick indices as age and metal-abundance 
indicators. 

The clusters of choice for this exercise are selected on the basis of availability of good 
quality age and metal abundance determinations (Table [3]). They span a wide range of metal- 
hcities, from [Fe/H] ~ -1.2 (NGC 6121) to nearly solar metallicity (M 67 and NGC 6528), 
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and include both old and intermediate-age (M 67) clusters. Lick indices for these clusters 
are presented in Table [21 



3.1. Reality Check 

In Figure S] cluster data are compared with the SOT solar abundance pattern models 
in a few index-index plots to illustrate that solar-scaled models cannot match the full set 
of observed cluster indices. The indices displayed are the default age and metal-abundance 
indicators employed by EZ_Ages. S07 showed that R(3, (Fe), Mg b, C24668, CNi (or CNa), 
and Ca4227 are the most reliable, recommended indices for age, [Fe/H], [Mg/Fe], [C/Fe], 
[N/Fe], and [Ca/Fe] determinations. This will be explored further in §3.2[ Henceforth, we 
will refer to this set of indices as the standard set. 

The models displayed in Figure S] were computed for the solar abundance pattern, and 
the data were measu red in the integrated spectra of the above gl obular clusters from th e 



Schiavon et al.l (120051 ) spectral library and that of M 67, obtained by lSchiavon et al.l (l2004al ). 
and revised by S07. In the upper left panel, model and data are compared in the (Fe) vs. H/? 
plane. S07 showed that these indices are excellent indicators of [Fe/H] and age, respectively. 
While the (Fe) average index is relatively insensitive to the abundances of elements other 
than Fe, H/3 is mostly sensitive to age, being only mildly sensitive to the overall metallicity 
of the stellar population. Therefore, the age and iron abundance of each cluster can be read 
fairly accurately from the position of the cluster data on the model grid — except for the case 
of NGC 6121 where the presence of blue HB stars affects the H/3-based age (see discussion 
in §3.2p . The remaining panels show plots of other metal-line indices against H/?, and (Fe) 
against B.6f. In the metal- line vs. H/? plots, one can see that the positions on the model 
grids where the data for most clusters fall are different from those they occupy on the (Fe) 
vs. H/9 plot. As discussed in §2.3[ this is a clear indication that the abundance pattern of 
the cluster stars is different from that of the Sun. We apply the method presented in the 
previous section to determine the abundance pattern of the clusters. 

We ran EZ_Ages on the cluster data presented in Table [2l The results are listed in 
Tabled together with abundance and age determinations from the literature. The literature 
values are based on abundance analyses of high-resolution spectra of individual cluster stars 
and on deep cluster CMDs, respectively. Differences between our results and those from the 
literature are displayed in Figure O In that plot, the dashed lines indicate the ±0.1 dex 
goal stated in the EZ_Ages specifications discussed in ^ Average values from the literature 
are displayed in black and individual determinations are shown in gray. For carbon and 
nitrogen, the values displayed are averages between CN-strong and CN-weak stars, as this 
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best approximates the effect of unresolved spectroscopic observations. 

A first inspection of tfie results shows that most abundance points fall between the two 
dashed lines, indicating that we are meeting the 0.1-dex goal for most of the abundances 
and for most clusters. Agreement is in fact particularly good for the abundances of carbon 
and nitrogen, where the 0.1-dex goal was reached for all clusters for which these elemental 
abundances are available in the literature. For magnesium and iron, we meet the 0.1-dex 
goal for some clusters, and the 0.2-dex requirement for all of them, with one exception: for 
NGC 6528, we fin d [Fe/H] = —0 . 26, in reas onably good agreem ent with the value [Fe/H] = 



—0.15 reported by lZoccali et al.l (120041 ) and [Bar 



3UV et al.l (120041) but in significant disagree- 



ment with the value [Fe/H] = +0.1 reported by ICarretta et al.l (1200 ll ). For calcium, there is 
one cluster only for which the 0.2-dex requirement is not met. For all the others, we meet 
the 0.1-dex goal. Given that this is the first systematic attempt at deriving abundances for 
all these five elements from integrated light, these results can be considered very satisfactory. 
Below we discuss the few exceptions for which agreement is not so good. 

Although the values of [Fe/H] determined by EZ_Ages are almost all within the 0.2-dex 
requirement, they are consistently slightly lower than the literature abundances determined 
from invidual stars. The zeropoint for [Fe/H] in the models is based upon the [Fe/H] mea- 
surements for the stellar library from which the S07 model is constructed (the S07 models 
have not been separately calibrated to the Solar spectrum) and the zeropoint is not known to 
perfect accuracy. A systematic 0.1-dex underprediction of [Fe/H] by the models is possible. 
However, this is not necessarily more uncertain than the zeropoints in the cluster abundance 
determinations from individual stars. In the case of NGC 6528, for which multiple litera- 
ture sources are available, the abundance determinations in the literature show substantial 
scatter: 0.25 dex in [Fe/H], 0.07 dex in [Mg/Fe], and 0.63 in [Ca/Fe]. In this context, it 
is difficult to address zeropoint issues on the order of 0.1 dex by comparison with existing 
cluster data. 

In general we note that, for some clusters, elemental abundances determined by different 
groups are in substantial disagreement, which makes the task of interpreting our results a 
little more complicated. That is especially the case of NGC 6528 mentioned above, for 
which iron abundances from two groups differ by ~ 0.25 dex, and calcium abundances differ 
by more than 0.6 dex! While a definitive statement on our ability to match that cluster's 
abundances with EZ_Ages awaits the solution of this discrepancy, we call attention to the 
fact that agreement with the average values is very good. 

EZ_Ages also performs very well in age determinations. For all clusters, with the sole 
exception of NGC 6121, the 0.2-dex requirement is met, and for 47 Tuc and M 67, the 
0.1-dex goal is also achieved. For NGC 6121, the age obtained is too young because of the 
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contribution by blue HB stars which tend to boost the strength of Balmer hnes in the spectra 



of old stellar populations, mimicking a younger age (see ^ 
on Balmer lines has been pointed out in previous studies ( 


3.21). The effect of blue HB stars 


Freitas Pacheco & Barbuv 
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2005: ISchiavon 


2007) and it 



can potentially lead to age underestimates. 



The cluster NGC 6441 also contains blue HB stars (IRich et al.l 119971 ) . but EZ_Ages 
determines a relatively old age for the cluster when using the H/3 index in the modelling 
process, which seems to contradict this interpretation. Unfortunately, we could not find 
a CMD-based age for this cluster in the literature so we cannot make a statement about 
how the age result from EZ_Ages compares to an independent age estimate. However, when 
the bluer Balmer line B.6f is used in the modelling of NGC 6441, EZ_Ages finds a younger 
age than that determined from H /j, suggesting that the i nfluence of the blue HB stars does 
show up in RSp. As discussed by lSchiavon et al.l (l2004bl ). the signature of blue HB stars is 
stronger at bluer wavelengths, so that they tend to impact age determinations based on H(5 
more strongly than those based on H/?. This can be clearly seen in the bottom right panel 
of Figure HI where both NGC 6121 and NGC 6441 look younger in the (Fe) vs. B.5f plot 
than in that involving H/3. 

In NGC 6441, the blue HB stars are greatly outnumbered by red HB stars, unlike 
NGC 6121 in which the blue HB population is substantial. They therefore have only a 
modest influence at the wavelength of H/3, but a larger influence at the blue r wavelength of 



B.6f, w here the blue HB stars contribute a greater fraction of the total light (ISchiavon et al. 
2004bl ). We further discuss this point in §3.2[ 



Finally, we note that the values we found for the abundances of C, N, and Ca for 
NGC 6121 are suspiciously low. Unfortunately, we could not find literature values for C and 
N abundances to compare with our estimates. However, the numbers Table [3] and Figure O 
show that our [Ca/Fe] estimate is too low by about 0.3 dex, which confirms our suspicions 
that there might be a problem with our procedure for [Fe/H] ^ —1.0. While this clearly 
deserves further investigation, for the time being we strongly discourage users from applying 
EZ_Ages in the [Fe/H] ^ —1.0 regime. 

Another check on the accuracy of the stellar population modeling is to compare the 
index predictions of the best-fitting model to the observed line strengths. Figure [6] shows 
the observed values of the Lick indices as measured in the five test cluster spectra, plotted 
against the predicted index line strengths from the best-fitting S07 model, as determined by 
EZ_Ages. Black and gray symbols show indices that are and are not included in the fitting 
process, respectively. The indices used in the fitting process are extremely well reproduced 
by the best-fitting model. In addition, the indices that are not used in the fitting process 
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are also reasonably well matched by the model predictions, with some scatter but with no 
indication of systematic problems in the modeling of individual indices. 

3.2. Consistency Check 

In this test, we verify how the results vary when different index sets are adopted as inputs 
to EZ_Ages. This is important, because the full set of indices considered by EZ_Ages are not 
always available to observers, due to limitations such as those determined by the instrumental 
setup adopted, or by the redshift of the sample studied. Therefore, it is important to 
make sure that age and metal abundances obtained from EZ_Ages do not depend on which 
absorption-line indices are employed — or, if they do, that the systematics is understood 
and can be accounted for. To test for this, we ran EZ_Ages on multiple combinations of 
line indices as measured for the clusters NGC 6441 and M 67, and inter-compared the 
ages and metal abundances obtained. As a basis for the comparisons, we adopt the results 
obtained when using the standard set, i.e., H/3, (Fe), Mg b, C24668, CNi, and Ca4227. The 
indices used in each separate model fit are listed in Table HI We substitute each alternative 
index separately to asses the effect due to that individual index. We also examine some 
combinations of indices, specifically a "balmer" model, in which we fit using the average of 
H/5, B.'jp, and B.5f, an "all" model, in which all lines are included in the fit, and a "high-z" 
model, which simulates the effect of applying EZ_Ages to higher redshift data, where only the 
bluer indices (A < 4400 A) may be available because the indices at longer wavelengths are 
redshift out of the optical spectrum. Results from the fitting process are provided in Table [5] 
and displayed in Figures [7] and [H where residuals relative to the age and metal abundances 
obtained with the standard index set are plotted for the various index combinations. 

Looking first at NGC 6441, Figure [7] shows that, for the vast majority of the index com- 
binations, ages and metal abundances vary by less than ±0.1 dex, indicating that EZ_Ages 
and the S07 models have attained an outstanding degree of consistency. There are a few ex- 
ceptions, though, which occur for the following index combinations: 1) When G4300 is used 
as a carbon abundance indicator in place of C24668 (panel c, triangles); 2) When higher or- 
der Balmer lines are used as age indicators (panel a); 3) When higher order Balmer lines are 
used as age indicators and Fe4383 is used as the iron- abundance indicator as in the "high-z" 
model (panel d, triangles) or when the higher order Balmer lines and Fe4383 are averaged 
in with the other lines as in the "all" model (panel d, squares). In case 1, the G4300-based 
[C/Fe] and [N/Fe] values differ from those obtained from the standard indices by -0.14 and 
+0.18, respectively. In case 2, ages tend to be younger when higher-oder Balmer lines are 
used in place of H/5, by up to ~ -0.2 dex. Finally, case 3 is similar to case 2, except that 
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[Fe/H] is too high, and consequently all the other abundances are too low. 

The difference found in case 1 when replacing C24668 by the G4300 index is due to the 
different sensitivities of the two indices to the ratio between the abundances of carbon and 
oxygen (C/0). The two indices are sensitive to variations of the oxygen abundance because 
of details of the molecular dissociation equilibrium in the atmospheres of cool stars. Of 
all molecules in whose formation carbon and oxygen take part, carbon monoxide (CO) is 
the hardest one to break, because it has the highest dissociation energy. Therefore, at the 
temperatures prevalent in the atmospheres of G and K stars, most available free carbon and 
oxygen atoms are locked in CO. As a result, variations in the abundance of oxygen, usually 
the most abundant of the two species, have a strong influence on the amount of carbon that 
is free to form other molecules, such as C2 and CH, which are responsible for the vibrational 
bands measured by the C24668 and G4300 indices, respectively. Therefore, an increase in the 
abundance of oxygen tends to provoke a decrease in the strength of those molecular bands. 

While the concentration of the C2 molecule depends quadratically on the abundance 
of carbon, that of CH depends only linearly on that abundance. As a result, C24668 is far 
more sensitive to carbon than to oxygen, while the the sensitivity of the G4300 index to 
those two abundances is approximately the same (though for both indices, the sensitivity to 
carbon and oxygen abundances have opposite signs). Although this makes the G4300 index 
a more uncertain indicator of carbon abundance, it raises the possibility that a combination 
of the two indices may be used to constrain the abundance of oxygen. We verified that 
possibility by raising the input oxygen abundance in the models from [0/Fe] = to +0.3. 
As a result, the carbon and nitrogen abundances inferred from use of C24668 and G4300 
indices agreed to within 0.05 dex. Furthermore, the ages and abundances inferred for the 
0-enhanced models were not substantially different from the results using the [0/Fe] = 
modelfl: [N/Fe] increased by +0.12 dex, while the ages inferred from and }i6p increased 
by 0.6 Gyr. The values of H/3 age, [Fe/H], [Mg/Fe], and [Ca/Fe] showed negligible change 
(< 0.1 Gyr or < 0.05 dex). Interestingly, the ox ygen abundance in NGC 6441 stars ranges 



from [0/Fe] = -0.05 to 0.34 (iGratton et al.ll2006l ). with stars in the upper end of the interval 



having "normal" oxygen abundances, and likely being more numerous and dominating the 
cluster light. Therefore, the value we obtained from the cluster's integrated light is in 



■^The changes in age and abundance determinations from increasing [0/Fe] to +0.3 dex are small because 
only the effect on the atmospheric line strengths are taken into account. The effect of enhancing [0/Fe] in 
the isochrone are not included because reliable 0-enhanced isochrones have not yet been incorporated into 
the S07 model. An 0-enhanced isochrone would result in younger age measurements without significantly 
changing the inferred abundances. S07 shows that model grids with a-enhanced isochrones move along 
lines of constant [Fe/H] in index-index plots. The result is that the inferred ages are different when using 
a-enhanced isochrones, but that the abundance determinations are relatively unchanged. 
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good agreement with the known cluster abundance, which makes it possible that oxygen 
abundances might be inferred from a combination of the C24668 and G4300 indices. This 
will be further investigated in a forthcoming paper. 

The dependence in case 2 of the resulting age on the Balmer line adopted is not sur 



prising, given that NGC 6441 is characterized by the presence of blue HB stars (IRich et al. 



1997 



Busso et al.ll2007l ). Because these stars have early-F and A spectral types, they tend to 
strengthen Balmer lines in the integrated spectra. That effect is obviously more important 
at bluer wavelengths, where the contribution by early-F and A stars to the integrated light 
is the greatest. Therefore, higher-order Balmer lines tend to be more strongly affected. As 
a result, when models that do not take into account blue HB stars are compared with the 
data, they tend to infer systematically younger ages, a nd mor e so according to higher order 



Balmer lines (see discussion in S07). ISchiavon et al.l (l2004bl ) have in fact shown that this 



effect can be used to constrain the morphology of the HB of globular clusters, solely on the 
basis of Balmer line strengths in their integrated spectra. In particular, the ratio between 
H5i? and H/5 was shown to be very sensitive to the presence of blue HB stars. 

Before discussing case 3, where Fe4383 and 116f are used in the fitting process, let us 
look at the results of the consistency test for M 67. Here again there is agreement to within 
±0.1 dex for any combination of indices used, with the exception of combining Fe4383 and 
HSp- The three Balmer lines all give consistent results for M 67 (panel a), because M 67 
does not contain blue HB stars and because the integrated spectru m of M 67 used in th e 



this analysis has been carefully constructed to exclude blue stragglers (ISchiavon et al.ll2004al ). 



Unfortunately, there is a problem (as y et not well-understood) with the G4300 measurements 



for this cluster (ISchiavon et al.ll2004al ). thus we cannot use M 67 to test our hypothesis about 



the effect of oxygen abundance on G4300 and C24668. 

Indeed, the only inconsistency in the abundance determinations for M 67 is in the "high- 
z" model, in which Fe4383 and }i6F are used in the fitting process (case 3 from above). Here, 
as with NGC 6441, we not only derive younger ages, but the resulting iron abundance is 
higher by ~ +0.24 dex. Consequently, all the other elemental abundances are found to 
be somewhat low. For M 67, we cannot invoke blue HB stars to explain this discrepancy. 
Interestingly, substituting B.6p or Fe4383 into the standard index set individually (triangles 
in panel a and squares in panel b) does not result in large discrepancies from the standard 
model. How then can we understand the dramatic differences seen when both H^^ and 
Fe4383 are used in the modelling process? 

Figure M shows solar abundance model grids from the S07 model. The M 67 data 
are overplotted as diamonds. In panel a, we see that H/3 and (Fe) produce model grids 
where lines of constant age and constant [Fe/H] are nearly perpendicular. This is one of 
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the reasons this index combination is so useful for age and [Fe/H] determinations. When 
Fe4383 is substituted for (Fe) (panel b), the lines of constant [Fe/H] are less nearly vertical, 
lessening the diagnostic power of the index-index grid. Likewise, when B.6f is substituted 
for H/5 (panel c), the lines of constant age are less nearly horizontal, again constraining the 
age and [Fe/H] determinations. This is because, in the blue region, line crowding is higher, 
making B.6f more sensitive to metallicity. Moreover, the blue spectral region is more strongly 
affected by warm stars than the red, making Fe4383 more sensitive to age than its redder 
counterparts. Notice, however, that there is still enough spread in the grids in panels b and 
c that small errors or zero-point uncertainties in the data or the models do not result in 
substatially different age and [Fe/H] measurements for M 67. 

Now consider panel d. When B.6f and Fe4383 are combined in the fitting process, the 
model grid lines are very far from perpendicular and the model space collapses down toward 
the familiar age-metallicity degeneracy. Here, small errors or zero-point uncertainties result 
in substantially different age and [Fe/H] measurements. Even the slight metallicity difference 
seen between panels a and b becomes a substantial discrepancy in panel d. This raises 
warning flags for using these indices to determine age and [Fe/H]. In particular, we note that 
Fe4383 appears to always give slightly higher values of [Fe/H] in the data presented both here 
and in ongoing work by the authors, suggesting that there is some inaccuracy in the modelling 
of this index. The fact that the discrepancy appears in the case of M 67, whose abundances 
are known to be solar-scale, suggests that this may be a zeropoint problem in the models, 
rather than a problem in calculating the abundance sensitivities of this index. Combined 
with the increased degeneracy of the H(5F-Fe4383 index-index space, this means that EZ_Ages 
should be used cautiously when determining stellar population parameters in higher redshift 
objects where H/3 is not available, and that the need for high S/N data is particularly acute 
in this regime. Even when the problem with Fe4383 is solved, the degeneracy of H(5ir-Fe4383 
index-index space means that small errors or zeropoint uncertainties in data can result in 
errors in the determination of age and [Fe/H] from these indices. This likely applies to all 
stellar population models, not just those presented in this work and S07. 

In summary, we conclude that EZ_Ages is characterized by a remarkable degree of 
consistency, in that ages and metal abundances inferred from different index sets agree in 
the vast majority of cases to within ±0.1 dex. Two of the inconsistencies we found, namely, 
that different Balmer lines indicate different ages, and different carbon-abundance indicators 
indicate different carbon abundances, can in fact be potentially explored to extract even 
more information from the integrated spectra of stellar populations. While blue Balmer 
l ines were shown to be useful to indicate the presence of blue HB stars in globular clusters 



( ISchiavon et al.ll2004bl ). the carbon index discrepancy can potentially be used to constrain 



the abundance of oxygen. The latter effect will be further investigated in a forthcoming 
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paper. The third inconsistency is due to the fact that the combined H(5i7'-Fe4383 index-index 
space is substantially more degenerate than H/3-(Fe) space, making accurate age and [Fe/H] 
determinations fr om these indices difficu lt. This is a matter of concern i n integrated studies o f 
globular clusters (ISchiavon et al.ll2004bl ) and higher redshift galaxies (ISchiavon et al.ll2006l ). 
for which only these bluer lines may be available. 



4. Comparison with Abundance Modeling of Galaxies 

Having performed the reality and consistency checks above, we now turn to comparisons 
between EZ_Ages results an d those coming from application of oth er models in the literature. 



We focus on the models by iThomas. Maraston. fc Benderl (120031 . hereafter TMB03), which 



also make predictions for variable abundance patterns. The TMB03 models also attempt to 
match non-solar abundance patterns by fitting Lick indices in unresolved stellar population 
spectra. These include both solar-scaled and a-enhanced models, in which the set of a 
elements N, O, Mg, Ca, Na, Ne, S, Si, and Ti are slightly enhanced, while the iron peak 
elements Cr, Mn, Fe, Co, Ni, Cu, and Zn are significantly depressed, in order to vary [a/Fe] 
at fixed total metallicity [Z/H]. In these models, [C/Fe] is fixed at solar. 



Thomas et al.l ( 120051 . hereafter T05) use the TMB03 models to estimate ages, [Z/H], 
and [a/Fe] for a sample of 124 early- type galaxies, based on measurement of the Lick indices 
H/3, (Fe), and Mg h. To compare the stellar population modeling results of using EZ_Ages 
and the S07 models to those using TMB03, we have run EZ_Ages on the sample of 124 
galaxies presented in T05, using the values of H/3, (Fe), and Mg h given in their Table 2. 
With these indices, EZ_Ages can determine the SSP age, [Fe/H], and [Mg/Fe] values of the 
sample galaxies and compare these results with those of the TMB03 models used in T05. 

Although the details of the model fitting procedure differ, T05 perform a very similar 
analysis to that of EZ_Ages, determining age, [Z/H], and [a/Fe] from the same Lick indices 
that are the standard set EZ_Ages uses to calculate age, [Fe/ H], and [Mg/Fe]. The models 



used in T05 use the same iKorn. Maraston. &: ThomasI (120051 ) index sensitivity functions as 



S07 and should therefore show the same index strength variations when the abundances are 
modified. One practical difference is that the TMB03 models are cast in terms of [Z/H], 
rather than [Fe/H]. In order to compare results from the two models, we convert T05's 
estimated [Z/H] into [Fe/H] using the conversion proposed by TMB03: [Fe/H] = [Z/H] 
-0.94 [a/Fe] (TMB03). 

The TMB03 models assume that all a-elements track Mg. EZ_Ages follows a similar 
convention, by default setting Na, Si, and Ti to track Mg in the modeling process, while Cr 
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follows Fe. A significant advantage of EZ_Ages and the S07 models is that they also track 
the [C/Fe], [N/Fe], and [Ca/Fe] abundance ratios using C, N, and Ca-sensitive Lick indices, 
rather than making assumptions about how these elements vary. However, these differences 
should not affect the comparison of age, [Fe/H], and [a/Fe] results presented here, as H/?, 
(Fe), and Mg b are predominantly sensitive to age, Fe and Mg only (see Tab l e [H), a nd are thus 



treated the same by both models through the iKorn. Maraston. fc ThomasI (120051 ) sensitivity 
functions. 

The results of fitting the T05 galaxies using EZ_Ages are shown in Figure [THl In each 
panel, the dashed line shows a one-to-one relation and the solid line shows the best fitting 
line to the data while keeping the one-to-one slope fixed. The resulting zeropoint offsets 
between the model fits are indicated in the upper left corner of each panel. In general, there 
is good agreement between the results coming from application of the two models to the 
same data set. The scatter is small and the zeropoint offsets are modest. The values of 
[a/Fe] agree to within a very small offset. The ages estimated in T05 are typically ~ 25% 
(0.13 dex) younger than those we find with EZ_Ages, while the T05 values of [Fe/H] are 
~ 0.08 dex higher. 

In addition to the zeropoint offsets in the age and [Fe/H] comparisons, the two different 
analyses do not truly follow a one-to-one slope for either parameter: the [Fe/H] relation is 
significantly steeper than the unity relation, while the age comparison shows some curvature. 
At old ages (log age ^ 1.0 as estimated by EZ_Ages) the slope is very nearly one-to-one, 
and the agreement between the two age estimates is much better, with a zeropoint offset 
of 15% (—0.07 dex) between the T05 ages and those from EZ_Ages. At intermediate ages 
(0.8 < log age < 1.0 as estimated by EZ_Ages) the relation is steeper than unity, while at 
young ages (log age ^ 0.8) the relation flattens out again, although at a substantial zeropoint 
offset. At intermediate and young ages, the T05 ages are 33% (0.17 dex in log age) younger. 
The origin of the curvature in the age relation is unclear. In any case, age agreement for 
old galaxies is very good (within 15%) and is still within 33% for all galaxies down to the 
youngest SSP ages estimated in the T05 sample. 

Based on the cluster comparison of §3.11 and Figure the zeropoint uncertainties in 
the S07 models are roughly 0.07 dex in log age, 0.1 dex in [Fe/H], and 0.05 dex in [Mg/Fe]. 
The differences between the T05 and the EZ_Ages results are therefore within the zeropoint 
uncertainties of the models, with the exception of the age estimates for galaxies younger 
than 10 Gyr. For these galaxies, the T05 results are younger by more than twice the indi- 
cated zeropoint uncertainties in the S07 and EZ_Ages models. The best way to resolve this 
discrepancy is to compare the models with data for metal-rich, intermediate-age clusters. 
While the models by S07 have been shown in the previous section and in Schiavon (2007) 
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to match the data for M 67 ( ~ 3.5 Gyr-old and solar metaUicity), more data are clearly 
needed to better constrain the models in this crucial age/metallicity regime. 

The non-unity slope of the [Fe/H] comparison also deserves further investigation. While 
the differences are negligibly small in the solar metallicity regime, they climb up to 0.2 
dex at the high-metallicity end. For the most metal-rich galaxies, EZ_Ages obtains [Fe/H] 
~ 0.2, and T05 find galaxies to be more metal-rich by ~ 0.2 dex. As in the case of age 
determinations, we again find ourselves lacking data that could help decide between the two 
model sets, given the non-existence of integrated spectroscopy for super-metal-rich clusters. 
It would be interesting, however, to compare results from the two models by confronting 
them with cluster data in a regime for which such data are available, such as moderately 
metal-poor clusters. 

In Figure [TT| we plot the difference between the values of [Fe/H] estimated in T05 and by 
EZ_Ages as a function of various other properties of the stellar populations . The top panel 
shows the [Fe/H] differences as a function of SSP age and a-enhancement. The differences 
in the [Fe/H] estimates do not appear to be strongly correlated with either SSP age or a- 
enhancement. In the lower panel, we show the differences in [Fe/H] as a function of the 
differences in the estimated ages and a-enhancements. The [Fe/H] differences are correlated 
with differences in the age estimates between the two different models. Furthermore, the 
correlation is in the direction expected from the correlated errors in the index-index plots: 
where T05 obtain substantially younger ages than EZ_Ages, they also find higher [Fe/H]. 
The difference in the age estimates by the different models therefore naturally explains the 
observed differences in [Fe/H], both the zeropoint shift and the non-unity slope of the [Fe/H] 
comparison. Estimates of [a/Fe] are substantially less biased by a zeropoint error in the 
age. That is because the Mg h and (Fe) indices have very similar age- dependence, so that 
the zero-point uncertainties in the parameters determined by these two indices ([Mg/H] and 
[Fe/H]) cancel out. 

From this analysis, it is clear that the TMB03 models and EZ_Ages give fairly consistent 
results, with a relatively small disagreement in the age zeropoints of the models which 
produces modest disagreements in [Fe/H]. Abundance ratios such as [a/Fe] are much less 
sensitive to zeropoint offsets. The comparison of EZ_Ages results with ages from CMD 
fitting of Galactic clusters indicates that the age zeropoint in EZ_Ages is correct to 
within 0.07 dex in log age. The [Z/H] and [a/Fe] results from TMB03 models have been 



compared with cluster data (jMaraston et al.ll2003l : iThomas. Maraston. &: Benderll2003l ). but 



age estimates have not been independently tested, therefore it is difficult to say what the 
zeropoint uncertainties in the age determinations for these models should be. We therefore 
conclude that EZ_Ages and the S07 models give results in good agreement with the TMB03 
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models (at least in the high-metallicity regime tested in this analysis of galaxy spectra), but 
that at intermediate ages (age < 10 Gyr) the TMB03 SSP ages may underestimate the true 
SSP ages by as much as 33% (0.17 dex in log age). 

A more general result from this comparison is that the sequential grid inversion algo- 
rithm presented in this paper, combined with the SSP models of SOT, does an excellent job of 
reproducing [Mg/Fe] abundances obtained by other groups with other models. Although the 
EZ_Ages abundance results for [C/Fe], [N/Fe], and [Ca/Fe] in unresolved stellar populations 
cannot be compared to other models (since no other models exist for the comparison), the 
modeling process for fitting these other abundances is the same as that used to fit [Mg/Fe] 
and the success of this process for [Mg/Fe] bodes well for other abundances. The test pre- 
sented in this section, combined with the stringent tests with cluster data described in ^ 
are reassuring evidence that this modeling process is reasonable and can be used to provide 
a quantitative assessment of multiple elemental abundances. 



5. Conclusions 



In this paper, we have presented a methodology for measuring SSP ages, [Fe/H], and 
individual abundance ratios [Mg/Fe], [C/Fe], [N/Fe], and [Ca/Fe] for unresolved stellar pop- 
ulations. We do this by exploiting the different sensitivities of various Lick indices to a 
variety of elemental abundances, and the ability of the S07 stellar population models to ac- 
curately model the Lick indices for a wide range of input abundance patterns. The algorithm 
presented here has been implemented in the IDL code package EZ_Ages, which is available 
for download and general use. 

We have subjected the modeling process described here to numerous rigorous tests and 
comparisons, with the following results: 



1. Comparison with Galactic cluster data: 

a. Ages: EZ_Ages age estimates reproduce the results of cluster CMD fitting to 
within 0.15 dex for all clusters, with the exception of NGC 6121, whose Balmer 
line strengths are affected by a blue HB. For two out of three of the remaining 
clusters, the EZ_Ages age is within 0.1 dex of the CMD result. 

b. [Fe/H] and [Mg/Fe]: EZ_Ages estimates match the results from high-resolution 
spectroscopy of individual cluster members to within 0.1 dex for some of the 
clusters, and to within 0.2 dex for all clusters (with the exception of the disagree- 



ment between our value of [Fe/H] for NGC 6528 and the results of ICarretta et al. 
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2001 



we are in ^ ±0.1 dex agreement with the results of IZoccah et aLll2004l and 



Barbuy et aLll2004j for [Fe/H] for this cluster). 

c. [C/Fe] and [N/Fe]: For the only two clusters where C and N abundances are 
available (47 Tuc and M 67), EZ_Ages results match those from high- resolution 
spectroscopy of individual cluster members to within 0.1 dex. For NGC 6121, we 
find suspiciously low values, that might indicate a problem at the low metallicity 
end of our models. Therefore, we strongly caution users against using EZ_Ages 
to determine these abundances for systems with [Fe/H] ^ —1.0. 

d. [Ca/FeJ: EZ_Ages results are within ~ 0.1 dex for all clusters except NGC 6121. 
The latter is certainly related to the problem mentioned above regarding [C/Fe] 
and [N/Fe] in the metal-poor regime, so we do not recommend use of EZ_Ages 
for [Ca/Fe] determinations for systems with [Fe/H] ^ —1.0. 

2. Consistency test performed on the cluster NGC 644 using different combinations of 
Lick indices in the fitting process yields results that are consistent to within ±0.1 dex 
for all stellar population parameters except for the following cases: 

a. Using G4300 instead of C24668 to fit [C/Fe] results in values of [C/Fe] and [N/Fe] 
that differ by —0.14 and +0.18 dex, respectively. This difference is probably due 
to differences in the index responses to O, which are not explored in this work. 

b. Using H7i? or ¥16 p instead of H/3 ge indicator results in younger ages for 
the cluster NGC 6441 due to the effect of the cluster's blue HB. This discrepancy 
does not exist for M 67, which does not have a blue HB, and whose integrated 
spectrum is constructed to exclude blue HB stars. 

c. Using the bluest Balmer and Fe lines available, HSp and Fe4383, can result in age 
and [Fe/H] measurements that are substantially different than those determined 
using the standard set of indices, due to the increased degeneracy of the model 
space for these indices. This poses a problem for high redshift studies of stellar 
populations, where often only the bluest lines are available for analysis. 

3. Error analysis test: EZ_Ages implements an algorithm that attempts to propagate 
errors through the modeling process in an efficient way that takes into account the 
dependence of abundance fitting on other parts of the population fitting process. Monte 
Carlo simulations show that this simplified error estimation does an excellent job of 
matching the true errors in the population parameters due to measurement errors in 
the Lick indices. 



4. Comparison with lThomas et al. hood ) results: 
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a. Results for [a/Fe] are consistent between EZ_Ages results and those using TMB03 
models. 

b. Age and [Fe/H] comparisons show httle scatter. There are zeropoint offsets of 
15%-33% (0.07-0.17 dex) in age and 0.08 dex in [Fe/H] between the two analyses, 
and a non- unity slope in the [Fe/H] comparison. The age differences are small for 
the oldest (age > 10 Gyr) galaxies and increase for younger galaxies. The [Fe/H] 
effects seem to be a natural result of the age differences. The S07 and EZ_Ages 
age zeropoints are shown here to be good to within 0.07 dex, while the TMB03 
age zeropoints are untested. This suggests that the ages derived in T05 may be 
too young by up to 33% (0.17 dex) for intermediate age galaxies (age < 10 Gyr). 

c. Absolute estimates of abundance ratios are more reliable than absolute estimates 
of age or [Fe/H] because zeropoint uncertainties tend to cancel out. 

d. The sequential grid inversion algorithm method is proven to work for estimating 
[Mg/Fe]. Since this method is the same as the one used to determine other 
abundance ratios, those estimates are likely to also be reliable, as demonstrated 
in the comparison with cluster data. 

Overall, EZ_Ages and the S07 SSP models do an excellent job of fitting the SSP age, 
[Fe/H], and abundance ratios [Mg/Fe], [C/Fe], [N/Fe], and [Ca/Fe] for unresolved stellar 
populations. The small zeropoint uncertainties in age and [Fe/H] estimates illustrated by 
comparisons with Galactic cluster data demonstrate that absolute estimates can be made for 
these quantities with high S/N data. Absolute estimates of elemental abundances are robust 
to the (small) zeropoint uncertainties. EZ_Ages and the S07 models therefore make it possi- 
ble, for the first time, to perform a quantitative assessment of multiple individual elemental 
abundances from medium-resolution spectra of unresolved stellar populations. Furthermore, 
the abundance fitting process can be run in an automated way on large data sets. With 
these innovations, stellar population analysis is better positioned than ever before to address 
the task of unravelling the star formation histories of stellar systems. 
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Table 1. Lick Indices in the S07 SSP Models. 



Index 


KMT Sensitivity* 


SWB Sensitivity* 


Balmer Indices 


HSf 
H7A 

H/3 


Fe, C" 
Fe, Mg+ 
C, Fe, Mg 
C, Fe 


Fe, C, V 
Fe, C, Si 
C, Ti, (Mg) 
C, (Si), (Mg) 
(Ni) 


Fe Indices 


Fe4383 
Fe5015 
Fe5270 
Fe5335 


Fe, Mg, C" 
Ti, Mg, Fe 

Fe 

Fe 


C, Fe 
C, Fe, (Ti) 
C, (Fe) 
Fe, C 


Mg Indices 


Mg2 

Mg b 


Mg, C 
Mg, Fe, C 


Mg, C, Fe 
Mg, Fe, (Cr) 


C, CH, and CN Indices 


CNi 
CN2 
G4300 
C24668 


C, N, 
C, N, 
C, 0, Fe+ 

c, 


C, N, 
C, N, 
C, Fe, Ti 
C, 0, (Si) 


Ca Indices 


Ca4227 


Ca, C 


Ca, 0, (CN) 


Note. 


— Index sensitivities as given by 



Korn. Maraston. fc Thomas jjo OS, KMT) and 
Serven. Worthev. fc BrilevI ||2005 . SWB). Only the 



top three element sensitivities are given, and only 
those with significance above la in the model 
spectra. For the KMT sensitivities, only those of 
turn-ofl^ and giant branch stars are shown, as these 
two components comprise 90% of the light in the 
S07 models. For the SWB sensivities, those under 
2cr are shown in parentheses. 

*In addition to the reported sensitivities for in- 
dividual element abundances, all indices listed here 
vary with total metallicity. 

"'"Indicates sensitivities that only appear at high 
metallicity. 
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Indicates sensitivities that only appear at low 
metallicity. 



Table 2. Selected Lick Indices for Clusters 



NGC 6121 47 Tuc NGC 6441 NGC 6528 M 67 GC Err* M 67 Err 



H<5f 


2.338 


0.682 


0.926 


0.298 


0.977 


0.117 


0.090 


H7f 


1.398 


-0.660 


-0.148 


-1.308 


-0.324 


0.215 


0.050 


H/3 


2.280 


1.602 


1.804 


1.572 


2.246 


0.083 


0.090 


Fe4383 


1.267 


2.491 


3.048 


4.666 


3.841 


0.267 


0.200 


Fe5270 


1.211 


1.913 


2.128 


2.799 


2.647 


0.142 


0.090 


Fe5335 


0.984 


1.685 


1.785 


2.417 


2.332 


0.079 


0.080 


Mg b 


1.577 


2.694 


2.651 


3.759 


2.927 


0.148 


0.100 


Mg2 


0.088 


0.159 


0.165 


0.249 


0.175 


0.015 


0.010 


G4300 


2.496 


4.656 


4.001 


4.819 


4.740 


0.198 


0.100 


C24668 


0.322 


1.680 


1.906 


4.510 


4.382 


0.131 


0.100 


CNi 


-0.070 


0.025 


0.020 


0.063 


-0.007 


0.010 


0.005 


CN2 


-0.047 


0.050 


0.047 


0.093 


0.018 


0.008 


0.008 


Ca4227 


0.321 


0.557 


0.581 


0.887 


0.927 


0.041 


0.100 



Note. — All indi ces are measured in A , except for CNi and CN2, which are measured in magni- 
tudes, as defined in lWorthev et ah I lll994h . 

'Uncertainties in the Lick index measurements for the Milky Way globular cluster (GC) data are 
dominated by the zeropoint uncertainties in converting measured indices to the Lick system and are 
thus the same for NGC 6121, 47 Tuc, NGC 6441, and NGC 6528. See ISchiavon et al.l for 
details. 
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Table 3. EZ_Ages Cluster Abundances Compared to Literature 







Age* 


[Fe/H] 


[Mg/Fe] 


[C/Fe] 


[N/Fc] 


[Ca/Fe] 


Ref. 


NGC 6121 


EZ_Ages 
literature 


' -"-0.7 

13* 


-1.30 ±0.09 
-1.18 ±0.00 


+0.36 ± 0.11 
+0.44 ± 0.02 


-0.12 ± 0.02 


0.15 ±0.15 


-0.05 + 0.10 
+0.26 + 0.02 


a, b 


47 Tuc 


EZ_Ages 
literature 


, „ „+maxt 
^•^■^-3.0 
12 


-0.80 ±0.09 
-0.7 ±0.05 


+0.28 ± 0.05 
+0.4 ±0.1 


-0.16 ± 0.03 
-0.2/0.0** 


+0.66 + 0.08 

+1.1/+0.3** 


+0.08 + 0.03 
+0.2 + 0.1 


c, d, e 


NGC 6441 


EZ_Ages 
literature 




-0.64 ±0.07 
-0.43 ±0.08 


+0.17 + 0.08 
+0.34 ± 0.09 


-0.20 ± 0.05 


+0.54 + 0.13 


+0.01 + 0.08 
+0.03 + 0.04 


f, g 


NGC 6528 


EZ_Ages 
literature 


14 g+maxf 
'-^■"-2.0 
11 


-0.26 ±0.06 
-0.15/+0.ltt 


+0.12 + 0.06 
+0.07/+0.14tt 


-0.04 ±0.04 


+0.31 + 0.08 


-0.07 + 0.05 
-0.407+0.23* 


h, i, j, k 


M 67 


EZ_Ages 
literature 


3.5 


-0.09 ±0.05 
0.0 ±0.1 


0.03 ±0.05 
0.0 ±0.1 


-0.07 ±0.03 
0.0 ±0.1 


-0.01 + 0.05 
0.0 + 0.1 


-0.03 + 0.07 
0.0 + 0.1 


1, m, n 



*A11 reported ages are estimated from H/3, except for 47 Tuc, which falls off the model grids in H/9-(Fe) and is therefore calculated 
using H7F. 

^Positive age errors cannot be calculated by EZ_Ages because they exceed the maximum 15.8 Gyr age of the models. 

'Age determined from relative ages of NGC 6121 and 47 Tuc in ISalaris fc Weis^ l l2002h . using 47 Tuc age from lSchiavon et al.l (l2002d) 

**CN-strong model / CN-weak model 

t*R esults differ substantially between authors. Reported abundances are lZoccali et all ll2004h and lBarbuv et al] ll2004h / ICarretta et al.] 
I2OOII 



References. — (a.) |Salaris & Weisdl2002l, (h) llvans et al.l 


19991. fc) ISchiavon et al.ll2002allbl (d) iBrilev et al.ll2004 (e) ICarretta et al.l 


2OO4I. ({). (s) iGratton et all ll2006l~). fhl lOrtolani et al.ll2001 


. (i) IZoccali et al.l 


2OO4I. fillBarbuv et alj|2004l. fk) ICarretta et ahlboOll. (11 


Schiavon et al.ll2004allbl. fm) iTautvaisiene et al.ll2000l. fn") IShetrone & Sanauist 


I2OOOI 



- 36 - 



Table 4. Models Using Different Indices For Fitting 



Indices used in fit 



Model name 


Balmer 


Fe 


Mg 


C 


N 


Ca 


Standard 


H/3 


(Fe> 


Mg b 


C24668 


CNi 


Ca4227 


H7F 




(Fe) 


Mg b 


C24668 


CNi 


Ca4227 


USp 


nSp 


(Fe) 


Mg b 


C24668 


CNi 


Ca4227 


B aimer 


all 


(Fe) 


Mg b 


C24668 


CNi 


Ca4227 


Fe4383 


H/3 


Fe4383 


Mg b 


C24668 


CNi 


Ca4227 


Fe5015 


H/3 


Fe5015 


Mg b 


C24668 


CNi 


Ca4227 


Fe5270 


H/3 


Fe5270 


Mg b 


C24668 


CNi 


Ca4227 


Fe5335 


H/3 


Fe5335 


Mg fe 


C24668 


CNi 


Ca4227 


Mg2 


H/3 


(Fe) 


Mg2 


C24668 


CNi 


Ca4227 


G4300 


H/3 


(Fe) 


Mg 6 


G4300 


CNi 


Ca4227 


CNa 


H/3 


(Fe) 


Mg fe 


C24668 


CN2 


Ca4227 


High-z 


USp 


Fe4383 




G4300 


CNi 


Ca4227 


All 


all 


all 


all 


all 


all 


Ca4227 



Note. — For indices used in the fitting process, bold face indicates the 
index that differs from the standard set. "All" indicates that all available 
indices for that element were averaged together (i.e., fitting for [C/Fe] was 
done using an average of C2 4668 and G4300). 



Table 5. EZ_Ages Results Using Different Indices For Fitting 



EZ_Ages results for NGC 6441 EZ_Agcs results for M 67 



Model 


Age 


[Fe/H] 


[Mg/Fe] 


[C/Fe] 


[N/Fe] 


[Ca/Fe] 


Age 


[Fe/H] 


[Mg/Fe] 


[C/Fe] 


[N/Fe] 


[Ca/Fc] 


Standard 


11.4 


-0.64 


+0.17 


-0.20 


+0.62 


+0.04 


4.0 


-0.09 


+0.03 


-0.07 


-0.01 


-0.03 


H7f 


6.9 


-0.54 


+0.17 


-0.22 


+0.60 


+0.05 


3.5 


-0.07 


+0.03 


-0.06 


-0.01 


+0.01 


nSp 


8.2 


-0.57 


+0.17 


-0.22 


+0.60 


+0.05 


3.7 


-0.08 


+0.05 


-0.06 


-0.01 


-0.04 


Balmer 


8.3 


-0.57 


+0.17 


-0.21 


+0.60 


+0.05 


3.7 


-0.08 


+0.03 


-0.06 


-0.01 


-0.01 


Fe4383 


11.1 


-0.55 


+0.08 


-0.25 


+0.58 


-0.04 


3.8 


-0.02 


-0.03 


-0.10 


-0.03 


-0.09 


FeSOlSt 














4.3 


-0.17 


+0.08 


-0.02 


-0.04 


+0.01 


Fe5270 


11.3 


-0.62 


+0.15 


-0.21 


+0.60 


+0.01 


3.9 


-0.10 


+0.03 


-0.07 


-0.01 


-0.03 


Fe5335 


11.5 


-0.66 


+0.19 


-0.19 


+0.64 


+0.06 


4.0 


-0.09 


+0.03 


-0.07 


-0.01 


-0.03 


Mg2 


11.4 


-0.63 


+0.22 


-0.20 


+0.60 


+0.02 


4.0 


-0.09 


+0.04 


-0.07 


-0.01 


-0.03 


G4300* 


11.6 


-0.65 


+0.15 


-0.33 


+0.79 


+0.02 














CN2 


11.4 


-0.64 


+0.16 


-0.20 


+0.56 


+0.02 


4.0 


-0.09 


+0.03 


-0.07 


-0.03 


-0.04 


High-z* 


5.7 


-0.34 


+0.03 


-0.33 


+0.50 


-0.10 


2.7 


+0.15 




-0.15 


-0.04 


-0.13 


All*t 


7.5 


-0.47 


+0.12 


-0.26 


+0.53 


-0.02 


3.7 


-0.07 


+0.03 


-0.06 


-0.02 


-0.01 



* There is a problem with G4300 in the M 67 spectrum jSchiavon et al ]|2004ah . thus for this cluster no separate model for G4300 is computed 
and the "high-z" and "all" models are computed with C24668. 

tFe5015 falls on a bad CCD column in the spectrum of NGC 6441, thus for this cluster no separate model for Fe5015 is computed and the 
"all" model does not include Fe5015. 
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1.0 1.5 2.0 2.5 3.0 3.5 1 2 3 4 
<Fe> Mg 6 




1.0 1.5 2.0 2.5 3.0 3.5 1 2 3 4 5 6 
<Fe> Mg 6 



Fig. 1. — Grids produced by the S07 model illustrating the effects of super-solar [Mg/Fe]. 
Solid lines show constant [Fe/H] from left to right of -1.3, -0.7, -0.4, 0.0, and +0.2 ([Fe/H] 
= 0.0 is shown as a dashed line for reference). Dotted lines show constant age from top to 
bottom of 1.2, 2.2, 3.5, 7.0, and 14.1 Gyr. The square shows an example data point and is 
the same in both panels. Upper panel: Models computed with solar-scaled abundances. The 
Balmer-Fe grid gives fiducial values oi t — 7 Gyr and [Fe/H] = -0.2. The Balmer-Mg b grid 
shows an [Fe/H] > 0.4 dex higher than the fiducial from the Balmer-Fe grid. This indicates 
that [Mg/Fe] is super-solar. Lower panel: Increasing [Mg/Fe] in the models slides the grid 
to the right toward the Balmer-Mg b data point, lowering the estimated [Fe/H]. B.f3 is also 
slightly affected by the abundance change, yielding a slightly younger fiducial age. In the 
Mg-enhanced model, the fiducial [Fe/H] from the Balmer-Fe grid matches the value of [Fe/H] 
estimated from the Balmer-Mg b grid, indicating that [Mg/Fe] = -1-0.3 is a good fit to the 
data. 
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Fig. 2. — Comparison of the error estimates determined by EZ_Ages with the results of 
the Monte Carlo error simulation. Age errors are plotted in the left panel and [Fe/H] and 
abundance errors in the right panel. The solid lines show one-to-one identity relations, 
with a ±30% spread indicated by dashed lines. The error estimates produced by EZ_Ages 
are consistent with those from the full Monte Carlo simulations, indicating that EZ_Ages 
produces reliable error estimates. 
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Fig. 3. — The distribution of EZ_Ages results for the Monte-Carlo error simulation. Yellow, 
red, blue, and green points correspond to the realizations of a single data point drawn from 
the 2%, 5%, 10%, and 20% error simulations, respectively. These trace out the error ellipses 
of the EZ_Ages fitting analysis and reveal the effect of correlated errors. 
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Fig. 4. — S07 solar-scaled model grids compared with cluster data. Model grids lines are as 
in Figure [H Because H/? and (Fe) are relatively insensitive to non-solar abundance ratios, 
the fiducial age and [Fe/H] values for each cluster can be read off the model grid in the H/3- 
(Fe) plot (upper left). If solar-scaled models were a good fit to all clusters, the data would 
fall in the same region of the model space in each index- index plot, as is approximately the 
case with M 67. The fact that other cluster data do not indicates that non-solar abundance 
patterns are needed to fit these clusters. 
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Fig. 5.— Comparison of EZ_Ages age (logt), [Fe/H], [Mg/Fe], [C/Fe], [N/Fe], and [Ca/Fe] 
estimates with values from the hterature. Error bars indicate the mean observational errors 
for our cluster data. Not all cluster parameters are available in the literature. For NGC 6528, 
two sets of abundances are available from different groups and these differe substantially (see 
Table [3]). The individual literature values are plotted in gray, with an average of the two 
literature values plotted in black. For 47 Tuc, [C/Fe] and [N/Fe] abundances in the literature 
are computed separately for CN-weak and CN-strong stars. Here we compare the EZ_Ages 
results for [C/Fe] and [N/Fe] against the average over both sets of stars, as this best simulates 
their combined contribution to the integrated light 
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Fig. 6. — Comparison of observed Lick index measurements with the predicted index values 
from the best-fitting S07 model, as determined by EZ_Ages. All indices included in the S07 
models are shown for each of the five test clusters. Measurement errors in the observed 
indices are approximately the size of the plotting symbols. Indices shown in black are those 
used in the abundance fitting process; those in gray are not used in the fitting. Sohd hues 
show the one-to-one relation. Indices not used in the fitting process show some discrepancies 
between model predictions and the observed values, but there is no evidence for systematic 
discrepancies in the modeling of any of the indices. 
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Fig. 7.— Variations in age (logt), [Fe/H], [Mg/Fe], [C/Fe], [N/Fe], and [Ca/Fe] estimates 
for NGC 6441 from using different indices in the fitting process. Model names correspond 
to tliose in Table HI Error bars indicate the observational errors calculated by EZ_Ages for 
NGC 6441 using the standard set of indices, (a-c) Results from substituting one alternative 
index into the standard index set (i.e., substituting G4300 for C24668 to fit [C/Fe]), and also 
for using averages of lines (i.e., an average of all Balmer lines instead of H/3). All combinations 
of indices give results consistent within ±0.1 dex, with the exception of the bluer Balmer 
lines (panel a), which underpredict the cluster age, and G4300 (panel c, triangle), which 
underpredicts [C/Fe] by —0.13 dex relative to the fit using C24668. The underprediction of 
[C/Fe] results in an overprediction of [N/Fe], as expected from the correlated errors between 
[C/Fe] and [N/Fe] (see Figure [3]). (d) Results for an average of all indices and for fitting 
when only indices blueward of 4400A are available (simulating higher redshift observations). 
Fitting with Fe4383 and H^^;' instead of (Fe) and H/? (as in the simulated high redshift case) 
results in substantially younger ages and higher [Fe/H]. See text §3.21 for details. 
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Fig. 8.— Variations in age (logt), [Fe/H], [Mg/Fe], [C/Fe], [N/Fe], and [Ca/Fe] estimates 
for M 67 from using different indices in the fitting process. Model names correspond to 
tliose in Table HI Error bars indicate the observational errors calculated by EZ_Ages for 
M 67 using the standard set of indices, (a-c) Results from substituting one alternative index 
into the standard index set (i.e., substituting G4300 for C24668 to fit [C/Fe]), and also for 
using averages of lines (i.e., an average of all Balmer lines instead of H/?). (d) Results for 
an average of all indices and for fitting when only indices blueward of 4700A are available 
(simulating higher redshift observations). Fitting with Fe4383 and instead of (Fe) and 
H/9 (as in the simulated high redshift case) results in substantially younger ages and higher 
[Fe/H]. See text gSjfor details. 




Fig. 9. — Solar abundance model grids showing the difficulty of measuring age and [Fe/H] 
from Fe4383 and B.6f. Lines of constant age run nearly horizontally and show (from top 
to bottom) 2.5, 2.8, 3.5, 5.0, 7.0, 8.9, and 14.1 Gyr models. Lines of constant [Fe/H] run 
nearly vertically and show (from left to right) —0.4, 0.0, and +0.2 dex models. The 3.5 
Gyr and solar abundance model are shown with dashed lines for reference. Measured index 
values for M 67 are overplotted as the diamond. The standard indices (panel a) produce 
model grids with lines of constant age and constant [Fe/H] nearly perpendicular. Line of 
constant [Fe/H] are more sloped from the vertical when using Fe4383 instead of (Fe) (panel 
b). Similarly, lines of constant age are more sloped from the horizontal when using B.6f 
instead of H/3 (panel c). Substituting one or the other of these indices for the standard set 
causes only small offsets from the standard results (see Figure [H]). However, the combination 
of Fe4383 and }16f (panel d) results in model grids that collapse down on one another, 
making the results much more sensitive to small errors or zeropoint uncertainties in either 
the models or the data. Because of the near-degeneracy of the models in Fe4383-H(5F space, 
the combination of these two indices can give results substantially different from the age and 
[Fe/H] measured using the standard index set. 
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Fig. 10. — Comparison of log Age, [Fe/H], and [a/Fe] results from T05 and from EZ_Agcs. 
The galaxy sample is that of T05, and the stellar population parameters are from SSP fits 
to H/9, (Fe), and Mg b for each set of models. T05 models are cast in terms of [Z/H], which 
we compare to EZ_Ages values of [Fe/H] using the conversion given in TMB03 for their 
a-enhanced mixture: [Fe/H] = [Z/H] - 0.94 [a/Fe]. We compare their values of [a/Fe] to the 
EZ_Ages results for [Mg/Fe], since T05 use the Mg b fine as their a-enhancement indicator. 
Dashed fines show the one-to-one relation. Solid lines show the best fit one-to-one slope 
allowing for a zeropoint offset. The size and direction of the zeropoint offset are indicated 
in the top left of each panel. Error bars in the lower right corner of each panel indicate 
the median observational errors from T05. Results for [a/Fe] are very consistent between 
the two models. T05 find ages which are younger by ~ 35% and [Fe/H] values which are 
higher by ~ 0.08 dex than the EZ_Ages results. Also, the slope of the [Fc/H] comparison 
differs somewhat from the one-to-one relation. These zeropoint differences, as well as the 
slope difference in [Fe/H], can be explained by differences in the model ages, which affect 
the [Fe/H] estimates, but have relatively fittle effect on [a/Fe]. The age zeropoint difference 
between the models is somewhat larger than the zeropoint uncertainty in the S07 models (the 
age zeropoint of the models used in T05 is uncafibrated) , suggesting that the T05 analysis 
may slightly underestimate the age of the sample galaxies. See text for details. 
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Fig. 11. — Top: Difference between [Fe/H] estimates for tlie T05 galaxy sample from EZ_Ages 
and from T05, as a function of stellar population age and [Mg/Fe]. Bottom: Difference in 
[Fe/H] as a function of differences in stellar population age and [a/Fe]. The differences 
between T05 and EZ_Ages estimates of [Fe/H] are most strongly correlated with differences 
in the age results between the two models (lower left panel). This suggests that the non- 
unity slope of the [Fe/H] comparison relation in Figure [10] is due to differences in the age 
estimates and arises from correlated errors (see discussion in §2.4p . The non-orthogonality 
of model grids in the index-index diagrams (as in Figure [1]) results in correlated errors such 
that underestimating the age of a population will result in overestimating the corresponding 
value of [Fe/H]. 



